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Abstract 

Recently in [1], the authors have given a method for calculation of the effective re- 
sistance (resistance distance) on distance-regular networks, where the calculation was 
based on stratification introduced in [2] and Stieltjes transform of the spectral distribu- 
tion (Stieltjes function) associated with the network. Also, in Ref. [1] it has been shown 
that the resistance distances between a node a and all nodes j3 belonging to the same 
stratum with respect to the a {Rap(i) , belonging to the i-th stratum with respect to 
the a) are the same. In this work, an algorithm for recursive calculation of the resistance 
distances in an arbitrary distance-regular resistor network is provided, where the deriva- 
tion of the algorithm is based on the Bose-Mesner algebra, stratification of the network, 
spectral techniques and Christoffel-Darboux identity. It is shown that the effective resis- 
tance on a distance-regular network is an strictly increasing function of the shortest path 
distance defined on the network. In the other words, the two-point resistance 
is strictly larger than The link between the resistance distance and random 

walks on distance-regular networks is discussed, where the average commute time (CT) 
and its square root (called Euclidean commute time (ECT)) as a distance are related to 
the effective resistance. Finally, for some important examples of finite distance- regular 
networks, the resistance distances are calculated. 

Keywords: resistance distance, association scheme, stratification, distance- 
regular networks, Christoffel-Darboux identity 

PACs Index: 01.55.+h, 02.10.Yn 
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1 Introduction 

A classic problem in electric circuit theory studied by numerous authors over many years is the 
computation of the resistance between two nodes in a resistor network (see, e.g., [3]). The study 
of electric networks was formulated by Kirchhoff [4] more than 150 years ago as an instance of 
a linear analysis. Besides being a central problem in electric circuit theory, the computation of 
resistances is also relevant to a wide range of problems ranging from random walks (see [5] ) , the 
theory of harmonic functions [6], to lattice Greens functions [7, 8, 9, 10, 11]. The connection 
with these problems originates from the fact that electrical potentials on a grid are governed 
by the same difference equations as those occurring in the other problems. For this reason, the 
resistance problem is often studied from the point of view of solving the difference equations, 
which is most conveniently carried out for infinite networks. In the case of Greens function 
approach, for example, past efforts [3], [12] have been focused mainly on infinite lattices. Little 
attention has been paid to finite networks, even though the latter are those occurring in real 
hfe. 

Within the theory of electrical networks, the standard method to compute the two-point 
resistances on a network is via the Moore-Penrose generalized inverse [13] or pseudo inverse 
of the Laplacian L of the underlying network, where the Laplacian is a matrix whose off- 
diagonal entries are the conductances connecting pairs of nodes. Just as in graph theory 
where everything about a graph is described by its adjacency matrix (whose elements is 1 if 
two vertices are connected and otherwise) , everything about an electric network is described 
by its Laplacian. It should be noticed that, the concept of two-point resistance called also 
effective resistance and resistance distance and molecular structure-descriptors based on them, 
were much studied in the chemical literature [14]- [25]. Also, it is closely related to the average 
first-passage time and the average commute time which are two important quantities in random 
walk models defined based on the Markov chain. It is shown in [26] that the computation of 
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the average commute time can be obtained from the pseudo inverse of the Laplacian L {L~^) 
of the underlying network. Also, it has been shown that this quantity and its square root is 
a distance, since it can be shown [26] that L"^ is symmetric and positive semidefinite. It is 
therefore called the Euchdean Commute Time (ECT) distance. In fact the ECT distance is 
the same as resistance distance or effective resistance (the effective resistance is symmetric and 
satisfies the triangle inequality and so is a distance metric). Therefore, any clustering algorithm 
(hierarchical clustering, /c-means, etc) which can be used in conjunction with the ECT distance, 
deals with the effective resistance. Also, Laplacian eigenmaps which is one of the graph- 
based approaches on dimensionality reduction and manifold learning and recently proposed by 
Belkin and Niyogi in [27] and share many similarities with other recent spectral algorithms for 
clustering and embedding of data, for instance Kernel PCA (Principle Component Analysis) 
algorithm [28] and spectral methods for image segmentation [29] (for a unifying view of the 
behavior of spectral embeddings and clustering algorithms, see [30]) deals with the Laplacian 
of the graph assigned to the data and uses the resistance distance. One of the most important 
aspects of spectral methods for clustering and embedding, including Laplacian eigenmaps, is 
the fact that they are all posed as eigenvalue problems. But clearly, for too large matrices, the 
computation by pseudo inverse becomes cumbersome. 

Recently, the authors have given a method for calculation of the resistance distance on 
distance-regular networks [1], where the calculation is based on stratification introduced in 
[2] and Stieltjes transform of the spectral distribution (Stieltjes function) associated with the 
network. Also, in Ref.[l] it has been shown that the resistances between a node a and all 
nodes P belonging to the same stratum with respect to the a {RajjH) , P belonging to the i-th 
stratum with respect to the a) are the same and the analytical formulas have been given for 
two-point resistances R^p(i),i — 1,2,3 in terms of the size of the network and corresponding 
intersection array without any need to know the spectrum of the pseudo inverse L~^. In this 
work, our starting point is along the same line by considering the Laplacian matrix associated 
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with a network, but we use the algebraic structure of distance-regular networks (Bose-Mesner 
algebra) such as stratification and spectral techniques specially the well known Christoffel- 
Darboux identity [31] from the theory of orthogonal polynomials to give a recursive formula 
for calculation of all of the resistance distances Ra/sii), i — 1,2, ...,d {d is the diameter of the 
graph) on the network without any need to the spectrum of the pseudo inverse L~^. In fact, 
in order to evaluate the resistance distance, one needs only to know the intersection array of 
the network. The other main result of the derived recursive formula is that, the resistance 
distance strictly increases by increasing the shortest path distance defined on the network, i.e., 
Raf3(m+i) — Rapim) > for all m = 1,2, ...,d — 1. 

The organization of the paper is as follows. In section 2, we give some preliminaries 
such as association schemes, distance-regular networks, stratification of these networks and 
Stieltjes function associated with the network. Section 3 is devoted to the concept of two- 
point resistance on regular networks and its relation with random walks. In Section 4 (main 
section of the paper) we show that the two-point resistance increases with increasing of the 
number of stratum and give a recursive formula for calculation of two-point resistances Ra/3(i) 
for i = 1, 2, ...d in distance-regular networks, where the algorithm is based on the Bose-Mesner 
algebra and spectral techniques specially the Christoffel-Darboux identity. Section 5 is devoted 
to calculation of effective resistances on some examples of distance-regular networks. The paper 
is ended with a brief conclusion and an appendix containing calculation of two-point resistances 
Rai3(*)j ^ = 1, 2, (i for some important finite distance-regular networks. 

2 Preliminaries 

In this section we give some preliminaries such as definitions related to association schemes, 
corresponding stratification, distance-regular networks and Stieltjes function associated with 
a distance-regular network. 
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2.1 Association schemes 

First we recall the definition of association schemes. The reader is referred to Ref.[32], for 
further information on association schemes. 

Definition 2.1 (Symmetric association schemes). Let F be a set of vertices, and let Ri{i — 
0, 1, d) be nonempty relations on V (i.e., subset ofVxV). Let the following conditions (1), 
(2), (3) and (4) be satisfied. Then, the relations {Ri}o<i<d onV xV satisfying the following 
conditions 

(1) {Ri}Q<i<d is a partition of V x V 

(2) Ro = {(«,«) :aeV} 

(3) Ri = R\ for < i < where R\ = {(/?, a) : {a, 13) E Ri} 

(4) For {a,j3) e Rk, the number p^j =\ {■j E V : (a, 7) G Ri and (7,/^) G Rj} \ does not 
depend on (a, /3) but only on i,j and k, 

define a symmetric association scheme of class d onV which is denoted hjY = (V, {Ri}o<i<d)- 
Furthermore, if we have p^j — p^^ for all i, j, /c = 0, 2, d, then Y is called commutative. 

The number v of the vertices, \V\, is called the order of the association scheme and R^ is 
called i-th relation (colour). For example, in the resistor networks the relations Ri, i = 0, 1, d 
can be interpreted as d + 1 different kinds of conductances, i.e., two nodes a, (5 have the i-th 
relation with each other if and only if the conductance between them be q (see Figure 1). In 
this paper we will deal with the special case where the conductance between two nodes a, /3 is 
c = ci if the nodes be adjacent, i.e., (a, /3) e Ri and the other conductances corresponding to 
the other relations will be taken to zero. 

Also note that, the intersection number p^^ can be interpreted as the number of vertices 
which have relation i and j with vertices a and /3, respectively provided that (a, pi) e Rk, and 
it is the same for all elements of relation R^. For all integers i (0 < i < d), set = and 
note that ACj ^ 0, since Ri is non-empty. We refer to Hi as the i-th valency of Y . 

For examples of association schemes, consider a cube known as Hamming scheme i?(3, 2), 
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in which V (the vertex set) is the set of 3-tuples with entries in F2 = {0, 1}. Two vertices 
are connected if and only if they differ by exactly one entry (see Figure 1). The distance 
between vertices, i.e. the length of the shortest edge path connecting them, will then indicate 
which relation they are contained in. E.g., if x — (0, 0,1), y — (0, 1, 1) and z = (1, 0, 1), then 
{x,y) e Ri, {x,z) e Ri and {y,z) G i?2- As an another example, consider the octahedron 
(a special case of a square dipyramid with equal edge lengths) which is the same as Johnson 
scheme J(4, 2), in which the vertex set V contains all 2-element subsets of the set {1, 2, 3, 4} 
and two vertices are adjacent if and only if they intersect in exactly one element. Two vertices 
are then at distance i, i = 0, 1, 2 if and only if they have exactly 2 — i elements in common 
(see Figure 2). 

Let Y — {X, {Ri}o<i<d) be a commutative symmetric association scheme of class d, then 
the matrices Aq, Ai, A^i defined by 



From (2-2), it is seen that the adjacency matrices Aq, Ai, Ad form a basis for a commutative 
algebra A known as the Bose-Mesner algebra of Y. Since the matrices A^ commute, they can 
be diagonalized simultaneously. The Bose-Mesner algebra has a second basis Eq, ...,Ed, such 
that, EiEj = 5ijEi and Ef=o -^j = ^ with Eo = 1/nJ (J is the all-one matrix) [32]. The 
matrices Ei, for (0 < i,j < d) are known as the primitive idempotents of the Y . Then, there 
are matrices P and Q such that the two bases of the Bose-Mesner algebras can be related to 
each other as follows 




(2-1) 



are adjacency matrices of Y such that 



d 



AiAj^Y.V%M- 



(2-2) 



d 



j=0 



1 " 



(2-3) 
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where, N denotes the cardinahty of X. Then clearly 



PQ = QP = NI. 



(2-4) 



It also follows that 




(2-5) 



which indicates that the Pjj is the i-th eigenvalue of Aj and that the columns of are the 
corresponding eigenvectors. Also, rrii :— trEi — N{a\Ei\a) is the rank of the idempotent Ei 
which gives the multiplicity of the eigenvalue P^j of Aj (provided that Pjj ^ P^j for k ^ i). 
Clearly, we have mo = 1 and J2f=Q rrii — N since, Z]f=o — ^- should be also noticed 
that, due to the tracelessness of the adjacency matrices Ai, i ^ 0, one can obtain the following 
identity (which we will use it in obtaining the main result of the paper) 



where, we have used the fact that {a\Ej\a) is independent of the choic of a & V. Finally the 
underlying network of an association scheme F = (V, Ri) is an undirected connected network, 
where the set V and Ri consist of its vertices and edges, respectively. Obviously replacing Ri 
with one of the other relations R^, i ^ 0, 1 will also gives us an underlying network F = {V,Ri) 
(not necessarily a connected network) with the same set of vertices but a new set of edges Ri. 

As we will see in subsection 2.3, in the case of distance-regular networks, the adjacency 
matrices Aj are polynomials of the adjacency matrix A, i.e., Aj — Pj{A), where Pj is a 
polynomial of degree j, then the eigenvalues Pij in (2-5) are polynomials of the eigenvalues 
Pji = Aj (eigenvalues of the adjacency matrix A). This indicates that, in distance-regular net- 
works (underlying networks of P-polynomial association schemes) the matrix P is a polynomial 



d d 



NSio = trAi ^N.Y, Pij{a\Ej\a) = ^ P/,m,- 



(2-6) 



j=0 j=0 
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transformation [33], as 



1 1 

Pi(Ao) Pi(Ai) 

P2(Ao) P2(Al) 



1 

■ Pi (Ad) 



\ 



(2-7) 



PdiXo) Pd{Xi) ... Pd{Xd) I 
or Py = Pj(Aj). For example, for the directed cyclic graph C„, where the adjacency matrices 
are given by Ai — S'^ with S as shift operator of period n {S'^ — /), we have Ai — and 
so the polynomial transformation in (2-7) reads as the well known Fourier transformation 



(the eigenvalues of the S are the n-th roots of unity, i.e., Aj 



a;' 



g27rj/n\ Pqj, undirected 



cychc graph C„, the adjacency matrices are given by A^ — S'^ + S~'^ and one can show that 
Ai = Ti{A/2), where Tj's are Chebyshev polynomials of the first kind. Then the polynomial 
transformation (2-7) is the same as the discrete cosine transformation. 



2.2 Stratifications 

For a given vertex a E V, let ri(a) :— {P E V : (a, /3) e Ri} denote the set of all vertices 
having the relation with a. Then, the vertex set V can be written as disjoint union of Ri{a) 
for i — 0,1, 2, d, i.e., 

V=(jUa). (2-8) 

1=0 

We fix a point o e as an origin of the underlying network of an association scheme, called 
reference vertex. Then, the relation (2-8) stratifies the network into a disjoint union of associate 
classes rj(o) (called the i-th stratum with respect to o). Let l'^{V) denote the Hilbert space 
of C-valued square-summable functions on V. With each associate class rj(o) we associate a 
unit vector in 1^{V) defined by 

\<t^i)-^ E (2-9) 
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where, \a) denotes the eigenket of a-th vertex at the associate class Ti{o) and Ki = |rj(o)| is 
called the i-th valency of the graph := p°j — \{'y : (0,7) e Ri}\ — |rj(o)|). The closed 
subspace of l'^{V) spanned by {|0j)} is denoted by A(r). Since becomes a complete 

orthonormal basis of A(G), we often write 

A(r) = ^©C|0,). (2-10) 

i 

Let Ai be the adjacency matrix of the network F = (V^, R). Then, from the definition of the 
i-th adjacency matrix Ai, for the reference state |0o) (|0o) = with o e as reference 
vertex), we have 

A|</'o)= E 1/5)- (2-11) 
Then by using (2-9) and (2-11), we have 

AM = ^/^^\(|>i). (2-12) 
2.3 Distance-regular graphs 

In this section we consider some set of important graphs called distance-regular graphs. First 
we recall the definition of so-called P-polynomial association schemes (which are closely related 
to the distance-regular graphs) as follows: 

Definition 2.2 (P-polynomial property) The symmetric association scheme Y — {X, {-Rj}o<i<(i) 
is said to be P-polynomial (with respect to the ordering Rq, ...,Rd of the associate classes) 
whenever for all i — 0,1, d, there exist di, e^, ff, di ^ ^ fi with: 

AiAi^diAi_i + eiAi + fiAi+i. (2-13) 

The condition (2-13) is similar to the well known three term recursion relations appearing 
in the theory of orthogonal polynomials, where Ai is in correspondence with x (see equation 
(2-29) in subsection 2.4). Using the recursion relations (2-13), one can show that A^ is a 
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polynomial in Ai of degree i, i.e., we have Ai = Pi{Ai) for i = 1,2, d. In particular, A = Ai 
multiplicatively generates the Bose-Mesner algebra (for more details see [34]). 

An undirected connected graph T = (y,Ri) is called distance- regular graph if it is the 
underlying graph of a P-polynomial association scheme, where the relations are based on 
distance function defined as follows: Let the distance between a,P eV denoted by d{a, (3) is 
the length of the shortest walk connecting a and (3 (recall that a finite sequence aQ,ai, ...,an G 
V is called a walk of length n if ctfe-i ~ for all k — 1,2, n, where ctfe-i ~ means that 
ctfe-i is adjacent with ak), then the relations Ri in distance-regular graphs are defined as: 
(a, P) e Ri if and only if d{a, /3) = i, for i = 0,1, ...,d, where d :=ma.x{d{a, (3) : a, (3 E V} is 
called the diameter of the graph. Since d{a, (3) gives the distance between vertices a and (3, d 
is called the distance function. Clearly, we have d{a, a) = for all a e F and d{a, /9) = 1 if 
and only ii a ^ j3. Therefore, distance-regular graphs become metric spaces with the distance 
function d. 

One should notice that, the condition (2-13) implies that for distance-regular graphs, we 
have the following relation 

Ti{(3) C r,_i(a) U r,(a) U r,+i(a), V /3 e r,(a). (2-14) 

We also note that, in distance-regular graphs, the stratification is reference vertex independent, 
namely one can choose every vertex as a reference one, while the stratification of more general 
graphs may be reference dependent. 

The relation (2-13) implies that in a distance- regular graph, p^^-^ = (for i Q, j dose not 
belong to {i — l,i,i + 1}), i.e., the non-zero intersection numbers of the graph are given by 

ai=Pli, bi=pl+i,i, Ci=p\_^^^ , (2-15) 

respectively (see Figure 3). The intersection numbers (2-15) and the valencies Ki satisfy the 
following obvious conditions 

+ + Cj = l^i l^i-lbi-l = KiCi, i = 1, ...,d. 
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Ko = ci = 1, bo = Ki = K, (co = bd = 0). (2-16) 

Thus all parameters of the graph can be obtained from the intersection array {bo, bd-i ; ci , . . . , c^} . 

By using the equations (2-2) and (2-16), for adjacency matrices of distance-regular graph 
r, we obtain 

AiAi = bi_iAi_i + QiAi + Q+iAj+i, i = l,2,...,d- 1, 

AiAd = bd-iAd-i + (« - Cd)Ad. (2-17) 
The recursion relations (2-17), imply that 



Ai = Pi{A), i = 0,l,...,d. 
By acting two sides of (2-17) on |0o) and using (2-12), we obtain 



(2-18) 



= + ^/Kiai\4)i) + i = 0, 1, d. (2-19) 

Then, by dividing the sides of (2-19) by and using (2-16), one can easily obtain the 
following three term recursion relations for the unit vectors |0i),i = O,l,...,d 



(2-20) 



where, the coefficients ctj and j3i are defined as 



ao = 0, ak = ak = K-bk- c^, ujf, = Pi = bk-iCk, k = l,...,d. 



(2-21) 



That is, in the basis of the unit vectors i = 0, 1, d}, the adjacency matrix A is projected 
to the following symmetric tridiagonal form: 

/ 



A^ 



ao Pi 
Pi Oil P2 
P2 eta P3 










... pd-i ad-i Pd 
... Pd ad J 



(2-22) 
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In Ref. [35], it has be shown that, the coefficients and Pi can be also obtained easily by 
using the Lanczos iteration algorithm. 

2.4 Stieltjes function associated with the network 

In this subsection we recall the definition of the Stieltjes function associated with a distance- 
regular network which is related to the spectral distribution corresponding to the network. To 
do so, first we recall some facts about the spectral distribution associated with the adjacency 
matrix of the network. In fact, the spectral analysis of operators is an important issue in 
quantum mechanics, operator theory and mathematical physics [36, 37]. Since the advent of 
random matrix theory (RMT), there has been considerable interest in the statistical analysis 
of spectra [38, 39, 40]. RMT can be viewed as a generalization of the classical probability 
calculus, where the concept of probability density distribution for a one-dimensional random 
variable is generalized onto an averaged spectral distribution of the ensemble of large, non- 
commuting random matrices. Such a structure exhibits several phenomena known in classical 
probability theory, including central limit theorems [41]. Also, the two-point resistance has a 
probabilistic interpretation based on classical random walker walking on the network. Indeed, 
the connection between random walks and electric networks has been recognized for some time 
(see e.g. [42, 43, 44] ), where one can estabhsh a connection between the electrical concepts of 
current and voltage and corresponding descriptive quantities of random walks regarded as finite 
state Markov chains (for more details see [5]). Also, by adapting the random- walk dynamics 
and mean-field theory it has been studied that [45] , how the growth of a conducting network, 
such as electrical or electronic circuits, interferes with the current flow through the underlying 
evolving graphs. In Ref. [35] it has been shown that, there is also connection between the 
mathematical techniques such as Hilbert space of the stratification and spectral techniques 
(which have been employed in [2, 46, 47, 48] for investigating continuous time quantum walk 
on graphs) , and electrical concept of resistance between two arbitrary nodes of regular networks 
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and so the same techniques can be used for calculating the resistance. Note that, although we 
take the spectral approach to define the Stieltjes function in terms of orthogonal polynomials 
(which are orthogonal with respect to the spectral distribution /i associated with the network) 
with three term recursion relations, in practice as it will be seen in the section 3, we will 
calculate two-point resistances without any need to evaluate the spectral distribution /i. 

It is well known that, for any pair [A, \4>o)) of a matrix A and a vector l^o), it can be 
assigned a measure /i as follows 

/x(x) = (0o|^(2:)|0o), (2-23) 

where E{x) — J2i \ui){ui\ is the operator of projection onto the eigenspace of A corresponding 
to eigenvalue x, i.e., 

A = j xE{x)dx. (2-24) 
It is easy to see that, for any polynomial P{A) we have 

P{A) = J P{x)E{x)dx, (2-25) 

where for discrete spectrum the above integrals are replaced by summation. Therefore, using 
the relations (2-23) and (2-25), the expectation value of powers of adjacency matrix A over 
starting site |0o) can be written as 

(0o|A"^|0o) = / x^i^idx), m = 0, 1, 2, .... (2-26) 

The existence of a spectral distribution satisfying (2-26) is a consequence of Hamburgers the- 
orem, see e.g., Shohat and Tamarkin [[49], Theorem 1.2]. 

Obviously relation (2-26) implies an isomorphism from the Hilbert space of the stratification 
onto the closed linear span of the orthogonal polynomials with respect to the measure /i. More 
clearly, the orthonormality of the unit vectors \4>i) implies that 

^ij = = -^(0o|A^,#o) = / P:{x)P^{x)f,{dx), (2-27) 



Effective resistances 



15 



where, we have used the equations (2-12) and (2-18) to write 

10.) = 4=^#o) = -^P^{AMo) = PKAMo), (2-28) 

with Pl{A) :— -^^Pi{A). Now, by substituting (2-28) in (2-20), we get three term recursion 
relations between polynomials Pj{A), which leads to the following three term recursion relations 
between polynomials Pj{x) 

xPl{x) = Pk+iPU{x) + akPl{x) + PkP'k-M (2-29) 

for A; = 0, 0? — 1, with Pq{x) ~ 1. Multiplying two sides of (2-29) by (3i...Pk we obtain 

P,...PkxPi{x) = P^...Pk+lPUli^) + akPi...PkPl{x) + PlA-Pk-iPLii^). (2-30) 

By rescaling Pj. as Qk — l3i...j3kPk, the spectral distribution /i under question is characterized 
by the property of orthonormal polynomials {Qk} defined recurrently by 

Qo{x)^l, Qi{x)^x, 

xQk{x) = Qk+i{x) + akQk{x) + (3lQk-i{x), k > 1. (2-31) 

It is a well known result from the theory of orthogonal polynomials that, the polynomials 
defined by (2-31) can be also evaluated via the following determinant 

x-ao 1 0...00 

LVi x-ai 1 0...00 

CU2 x-a2 I ... 



Qk{x) 



(2-32) 



... uJk-2 x-ak-2 1 

... Uk-i x-ak-i 

If such a spectral distribution is unique, the spectral distribution // is determined by the 
identity 

1 Q2i(:^) 



G,{x) 



H{dy) 



R X — y X — ao 



x—cci 



Qd{x) 



1=0 X Xi 



(2-33) 
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where, xi are the roots of the polynomial Qd{x). The function G^{x) is called the Stielt- 
jes/Hilbert transform of spectral distribution jj, or Stieltjes function and polynomials {Q^^"*} 
are defined recurrently as 

xQ^^\x) = g«,(x) + «,+ig«(x) + (5U,Q^il,{x), k > 1, (2-34) 

respectively. The coefficients bi appearing in (2-33) are calculated as 

bi^ \im{x-xi)G^{x). (2-35) 

Now if G^{x) is known, then the spectral distribution can be recovered from Gn{x) by 
means of the Stieltjes inversion formula: 

1 fy 

fj,{y) — fj,{x) — limy^o+ / ImGn{u + iv)du. (2-36) 

TT Jx 

Substituting the right hand side of (2-33) in (2-36), the spectral distribution can be determined 
in terms of x/, / = 1, 2, the roots of the polynomial Qd{x), and Guass quadrature constants 
bi, 1 — 1,2, ... as 

li{x) = Y.bi5{x-xi) (2-37) 
I 

(for more details see Refs.[49, 50, 31, 51]). Note that, by using (2-23) and (2-37), we have 

1 TTl ' 

h = lJi{xi) = {a\E{x,)\a) = -trE^ = (2-38) 
2.4.1 The ChristofFel-Darboux Identity 

In the next section we will treat with the calculation of the effective resistances in distance- 
regular networks, where we will use one of the most important theorems of the orthogonal 
polynomials known as Christoffel-Darboux Identity. This identity is expressed as follows 
Theorem (ChristofFel-Dcirboux Identity) Let {Qn{x)} satisfy (2-31). Then 

^Qk{x)Qk{u) , .-iQn+l{x)Qn{u) - Qn{x)Qn+l(u) 

UJiUJ2...UJk X-U 

For the proof, the reader is refered to [31]. 
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3 Two-point resistances in regular resistor networks 

A classic problem in electric circuit theory studied by numerous authors over many years, is 
the computation of the resistance between two nodes in a resistor network (see, e.g., [3]). The 
results obtained in this section show that, there is a close connection between the techniques 
introduced in section 2 such as Hilbert space of the stratification and electrical concept of resis- 
tance between two arbitrary nodes of regular networks and these techniques can be employed 
for calculating two-point resistances. 

For a given regular graph F with v vertices and adjacency matrix A, let r^j = r^j be the 
resistance of the resistor connecting vertices i and j. Hence, the conductance is Cij — — cji 
so that Cij = if there is no resistor connecting i and j. Denote the electric potential at the 
i-th vertex by and the net current flowing into the network at the i-th vertex by /j (which 
is zero if the i-th vertex is not connected to the external world). Since there exist no sinks 
or sources of current including the external world, we have the constraint J2^=i -^i = 0. The 
Kirchhoff law states 

Cij{Vi-V,) = I„ i = l,2,...,v. (3-40) 

Exphcitly, Eq.(3-40) reads 

LV = I, (3-41) 

— * — * 

where, V and / are T;-vectors whose components are and /j, respectively and 

L = ^c,|i)(i|-^c,,.|i)(j| (3-42) 



is the Laplacian of the graph F with 



Ci= Yl ^ij^ (3-43) 



for each vertex a. Hereafter, we will assume that all nonzero resistances are equal to 1, then 
the off-diagonal elements of —L are precisely those of A, i.e., 

L^kI-A, (3-44) 
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with K = deg{a), for each vertex a. It should be noticed that, L has eigenvector (1, 1, 1)* 
with eigenvalue 0. Therefore, L is not invertible and so we define the psudo-inverse of L as 

= Y: K'E^, (3-45) 

where, Ei is the operator of projection onto the eigenspace of corresponding to eigenvalue 
Aj. It has been shown that, the two-point resistances Rap are given by 

R^p = {a\L'^a) + {l3\L-^\p) - {a\L-^\l3) - {l3\L-^a). (3-46) 

This formula may be formally derived using Kirchoff 's laws, and seems to have been long known 
in the electrical engineering literature, with it appearing in several texts, such as Ref.[52]. For 
distance-regular graphs as resistor networks, the diagonal entries of L'^ are independent of 
the vertex, i.e., = L^j for all a,l3 & V. Therefore, from the relation (3-46) and the fact 
that is a real matrix, we can obtain the two-point resistance between two arbitrary nodes 
a and (3 as follows 

^a/5 = 2(L-^-L-^). (3-47) 

It should be noticed that, in distance-regular graphs (more generally, for underlying networks 
of association scheme) due to the stratification of the network, all of the nodes belonging to the 
same stratum with respect to the reference node, i.e., a, possess the same two-point resistance 
with respect to the a. More clearly, for all /3 e r^(Q;) we have 

L;^(^) = {a\L-'\P) = -L(a|L-Vm) = — (3-48) 

Then, by using (3-47), we obtain 

2 2 K 1 — P (A) 

^a/3W = — {«^m^aa ' (^m^"')aa} = — (a| . > ) , V /? G !„(«). (3-49) 

where, the upper index m in and Rapim) indicate that (3 belongs to the m-th stratum 

with respect to a. 
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As, Eq.(3-47) implies, in order to evaluate two-point resistances Rapim) , we need to calculate 
the matrix entries and L~p. To this end, one needs to know the spectrum of the pseudo 



inverse (see Eq.(3-45)) which is a task with high complexity for networks with large size, 
even with computer. In the following we give an algebraic method such that the two-point 
resistances are calculated recursively without any knowledge about the spectrum of the pseudo 
inverse of Laplacian of the network. 

3.1 Random walks and electrical networks 

The computation of effective resistances is relevant to a wide range of problems ranging from 
random walks (see Ref. [5]). Random walks on graphs are the bases of a number of classical 
algorithms. Examples include 2-SAT (satisfiability for certain types of Boolean formulas), 
graph connectivity, and finding satisfying assignments for Boolean formulas. It is this success 
of random walks that motivated the study of their quantum analogs in order to explore whether 
they might extend the set of quantum algorithms. In Refs. [2, 46, 47, 35, 48], the same 
techniques introduced in this paper in order to evaluate the effective resistances, such as 
the algebraic structure of distance- regular graphs (Bose-Mesner algebra), stratification and 
spectral analysis methods have been used for investigation of the continous time quantum 
walks on the regular networks. In order to show this connection more clearly, in the following 
we discuss the link between resistance distance and two important quantities (average first 
passage time and average commute time) defined in random walks on graphs. 

Let r be a complete undirected graph with N vertices numbered 1, 2, A^, in which each 
edge {a, /3) is assigned its distance d{a, /3) = da/j > 0. One can study the harmonic random 
walk in F with escape probability (the probability that a walk starting at a reaches /3 before 
it returns to a) 



In the literature, the harmonic random walk was often studied using techniques from elec- 




(3-50) 
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trical network theory. Denote by Ha/s, the hitting cost of the harmonic walk from a to (3, 
defined as the expected cost (total distance) to reach /3 for the first time when started from 
a. By elementary probabihty, the costs of the random walks ending at vertex N satisfy the 
following system of equations 



and Hn^n = 0. As noted in [5], the hitting costs have an interpretation in terms of electrical 
networks. We can think of F as an electrical network in which each edge {a, (3) has resistance 
dap (so, the Eq.(3-50) is a probabihstic interpretation of the effective conductance). If we inject 
current of value — 1 into each node and draw current of value N[N — 1) from node A^, then 
the voltages relative to node N estabhshed at the nodes satisfy the same equation as (3-51). 
Therefore the voltage at a is equal to the hitting cost Hap. Moreover, for a random walker on 
a network, one can assign a quantity known as average commute time (CT) denoted by n{a, j3) 
which is defined as the average number of steps the random walker, starting in state a ^ (5, 
will take before entering a given state (3 for the first time, and go back to a. Clearly, the 
average commute time is symmetric and is equal to n{a, (3) — m{f3\a) -\- m{a\f3) , where m{(3\a) 
is the average first-passage time defined as the average number of steps the random walker, 
starting in state a, will take to enter state (3 for the first time. Note that for random walks on 
distance-regular graphs, we have m{(3\a) — m{a\(3) and so, n{a,(3) — 2m{(3\a). By viewing 
the graph as an electrical network, the average commute time has an electrical equivalent 



where, N.k is the volume of the graph (the volume of a graph is defined as JZaev^a, with 
da as the degree of the vertex a). The Eq.(3-52) indicates that, the average commute time 
and effective resistance basically measure the same quantity. This quantity can also be called 
resistance distance (it has be shown that n{a, (3) is a distance measure). Further connections 
between random walks and electrical networks were explored by Chandra et al.[53]. In the 



Hai3 = XI Pocnidap + -f^/3,iv), for « 7^ N, 



(3-51) 



(3-52) 
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following we introduce a method for recursive calculation of resistance distances on distance- 
regular resistor networks based on spectral techniques specially by employing the Christoffel- 
Darboux identity. 

4 Recursive calculation of resistance distance based on 
spectral analysis methods and Christoffel-Darboux iden- 
tity 

In this section, we show that the resistance distance on distance-regular resistor networks 
increases with the number of the strata, i.e., is strictly larger than Ra/3(m) for m — 

1,2, ...,d — 1. In previous work [1], explicit formulas for the two-point resistances up to the 
third stratum, i.e., Ra/3(rn) for m = 1,2,3, have been given in terms of the intersection array 
of the network, where the authors have been employed the properties of the Stieltjes function 
associated with the network. Here in this work, we use the spectral techniques and Christoffel- 
Darboux identity and give a recursive formula for calculating the two-point resistances Rapim) 
ior m — 1,2, ...,d. In addition to the preference that this formula enables us to calculate 
recursively all of the resistance distances, it also indicates that the resistance distance strictly 
increases with the number of the stratum. 

Let a and (3 be two arbitrary nodes of the network such that /3 belongs to the m-th stratum 
with respect to a, i.e., (3 e T^i^a) (we choose one of the nodes, here a, as reference node). 
Now, for calculating the matrix entries and L^^ in (3-46), we use the spectral techniques 
to obtain 

^.--=h^w^|:h^w4i:^. (4-53) 

and 
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1 y^j ^^^ ^rnjK) Ej _ 1 y^miPm{\i) (4 54) 

where, we have considered Xq — k {k is the eigenvalue corresponding to the idempotent £^o)- 
Then, by using Eq.(3-47) and the fact that is independent of m (the number of stratum), 

we have 

2 ^ Km+lPm{^j) ~ l^mPm+l{Xj) _ "^y^ Hmf^m+l ^ Qm+l{l^)Qm{^j) ~ Qm{l^)Qm+l{^j 

NKm^rn+lj^^ ^ 1^ - Xj ^^«m'tm+l/3i-"/5m/5m+l ~^ 1^ - ^'j 

(4-55) 



Now, from the fact that {Qk{x)} satisfy the three-term recursion relations (2-31), we can 
the Christoffel-Darboux identity to write the right hand side of (4-55) as follows 

r.h.s — 2^ mj[u!i...u!m —) — 

r\ m d r\ m 

N ,^ ^ ^ E E ^.^KA,) = UNSio - ^i), (4-56) 

where, we have used the fact that 

= ^— il(fi;) = i^—f^h and 



and have done the following simplification by using (2-16) 

In the last equality of (4-56), we have used the distance- regularity of the network to substitute 
Pi{Xj) = Pij and then use the Eq.(2-6) (recall that mo = 1 and Pi{Xo) = Pi{i<-) = After 
these simplifications, we obtain the main result of the paper as follows 

2 ™ 

R^p(m+i) - R^p^m) ^ — — — (iV-^K;) > 0, m = l,2,...,d-l. (4-57) 

1\ KfYiOm i—Q 
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It should noticed that, by using the formula (4-57) we can evaluate the effective resistance 
between any two nodes recursively, if we know the two-point resistance RapW- In order to 
calculate Rapw, one can use the spectral techniues introduced in subsection 2.4 (equations 
(2-26), (2-37) and (2-38)) to write 



kI — A Jr-{k} k — X N .-^j, K — Xi 



and 

1^1 A I ^ If dii{x) 1 ^ rUiXi 



'kJ-A' ' kJr-{k}K-X NK.^^K-Xi ^ ' 



Then, by using (3-47), one can obtain 

2 „ mi{K-Xi) 2 „ 2(jV-l) 

- m S - ml .§„-"■ - ^v^- 

Note that, as the main result of the paper, the result (4-57) shows that the resistance 
distance (and consequently the other quantities related to the resistance distance such as the 
average first passage time and Euclidean commute time associated with a random walk) on 
distance-regular networks is strictly increasing function of the shortest path distance defined 
in subsection 2.3, i.e., the nodes belonging to the farthest stratum with respect to a, possess 
the smallest effective resistance with a. Apart from this fact, the formula (4-57) together with 
(4-60) gives an algebraic method for calculation of the resistance distances on distance-regular 
networks, where one needs only to know the intersection array of the networks without any 
knowledge about the spectrum of the pseudo inverse of Laplacian of the networks. 



5 Examples 

In this section, we calculate the effective resistances on the examples of Cycle network, d- 
cube network and Johnson network recursively by using the formula (4-57), where for the first 
example (cycle network), the general formula for the effective resistances Ra/3(i), i — 1, 2, d 
is deduced, whereas for the two latter ones the effective resistances are calculated up to the 
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third stratum. The effective resistances Rapi'): i = 1,2, ...,0? on some other important finite 
distance-regular networks is given in the appendix. 

5.1 Cycle network Cn 

The graph Cn ioi N = 2m or N = 2m + 1 consists of m + 1 strata. The intersection arrays 
for even and odd number of vertices is given by 

{bo, bm-i; ci, Cm} = {2, 1, 1, 1; 1, 1, 2} (5-61) 

and {6o,--,J'm-i;ci,...,c^} = {2,1,...,1;1,...,1,1}, (5-62) 

respectively. We consider the even case = 2m, the odd case can be considered similarly. 
For this case, we have kq — Km — 1; ki — 2, for Z = 1, ...,m — 1. Then by using (4-60) and 
(4-57), we obtain the effective resistances recursively as follows 

2m - 1 



2m - 3 _ 2(2m - 2) 
2m 2m 



2m -5 3(2m-3) , , 

i?a/J(3) = R.pi^) + = 2m ' - ^^"^^^ 

Prom (5-63), one can easily deduce the following result 

it:.,(0 = ^^^1^, l^l,2,...,m. (5-64) 

The above formula indicates that in the limit of the large number of vertices, i.e., in the limit 
m — > 00, where the cycle network tends to the infinite hne network, the effective resistances 
are given by 

Rapw^l: Z = l,2,... . (5-65) 
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5.2 d-cube 

The (i-cube, i.e. the hypercube of dimension d, also called Hamming cube, is a network with 
2*^ nodes, each of which can be labeled by an ci-bit binary string. Two nodes on the hypercube 
described by bitstrings x and y are are connected by an edge if \x — y\ — 1, where \x\ is the 
Hamming weight of x. In other words, if x and y differ by only a single bit flip, then the two 
corresponding nodes on the graph are connected. Thus, each of the 2*^ nodes on the d-cube 
has degree d. For the d-cube we have d + 1 strata with 

d' 

i^i^ ' ... , 0<i<d-l. (5-66) 
The intersection numbers are given by 

bi — d — i, < i < d — 1; Ci = i, 1 < i < d. (5-67) 
Then by using (4-60) and (4-57), we obtain 

2'^-1-d _ 2'^-i - 1 

Rapi^) - Ra^i^) + ^(^_ 1)2.-1 - (rf_l)2<i-2' 

2'^ -I- d- d{(l -l)/2 _ 1 2'\d^ -2d^2) -M{d-l) -2 

^a/3(3) - ^a/3(i)+ 2<^-''d{d - l){d - 2) ~ d{d - l){d - 2)^ 2^^ 

(5-68) 

The other resistance distances -Ra/jCi) for i > 4, can be calculated similarly. 

5.3 Johnson network 

Let n>2 and d < n/2. The Johnson network J(n, d) has all d-element subsets of {1, 2, n} 
such that two d-element subsets are adjacent if their intersection has size d—1. Two d-element 
subsets are then at distance i if and only if they have exactly d — i elements in common. The 
Johnson network J(n, d) has = d\{n^d)\ vertices, diameter d and the valency k — d{n — d). 
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Its intersection array is given by 

bi = {d-i){n-d-i), < i < d - 1; Ci = i'^, l<i<d, (5-69) 

Then by using (4-60) and (4-57), one can obtain 

_ 2{n\ - d\{n - d)\) 
~ din - d)n\ ' 

2[n\ - [l + d{n - d))d\{n - dy] _ 

c u ,^ r ^^ dUn - d)\(n - 2 - 2d(n - d)) . 
-{d{n -d)-{n-2) + '-^ ^ 



d{d — l){n — d){n — d — 1) n\ 

2 

R„p(3) = R^,3(2) H — — - — — — — — — — — - — - — — - — - — — — -{ 



d\{n-dy. 



d{d - l){d - 2){n - d){n -d- l)(n -d- 2)/4 'd!(n - d)\ 



d(n-d-d(d-l)(d-2)(.-d)(W-l)(W-2)/4)} ^ _2) X 
{d\n - 2d+ 1) + (3n - 2d{n - d) - iq) "^(^ ~ '^)'^;(^ ~ '^)' + _ _ _ ^wg^ _ 

4(d - l)(n -d-l) + 2(n- 2)(n - 4)1(1 - ^ , '' )}. (5-70) 
Again, one can obtain the other resistance distances Ra/3(i) for i > 4, similarly. 



6 Conclusion 

Based on the Bose-Mesner algebra corresponding to distance-regular networks, stratification, 
spectral techniques and Christoffel-Darboux identity, a recursive formula for calculating resis- 
tance distance in distance-regular resistor networks was obtained such that one can evaluate 
the resistance distances on these networks only by knowing the corresponding intersection ar- 
ray , without any need to know the spectrum of the pseudo inverse of the Laplacian of the 
networks. As an important result, it was shown that the resistance distance on a distance- 
regular network is an increasing function of the shortest path distance defined on the network. 
Although we focused specifically on distance-regular networks, we hope that the introduced 
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method might then be apphed to other underlying networks of association schemes which are 
not distance-regular ones such as finite and infinite square lattice and underlying networks 
of the root lattices of type An particularly finite and infinite hexagonal networks (n — 2) in- 
troduced in [47] by employing the Krylov-subspace Lanczos algorithm [35] iteratively to give 
three-term recursion relations to the networks, where these problems are under investigation. 
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Appendix 

In this appendix, we give the two-point resistances Rapii), i — 1,2, for some important 
finite distance-regular networks. 

1. CoUineeirity graph, gen. octagon (s, 1), s = 2,3,4 



Then, 



TV = / + 2s-^ + 2s^ + 2s + 1, {6o, &2, h] Ci, Ca, C3, C4} = {2s, s, s, s; 1, 1, 1, 2}, 

K — Oq — 2s, K,2 = = 2S , K3 — — 2s , K,4— — S . 

C2 C3 C4 



_ + 2s^ + 2S + 2 
- 54 + 2s3 + 2s2 + 2s + l' 
_ + 2s + 2 _ s3 + 3s2 + 4s + 4 

^a/3(2) - Ra/3W + ^4 + + 2s2 + 2s + 1 ~ S^ + 2s3 + 2s2 + 2s + 1 ' 

sJ-2 s^ + 3s^ + 5s + 6 

^a/3(3) - ^a/3(2) + ^4 + 2s3 + 2s2 + 2s + 1 ~ S^ + 2s3 + 2s2 + 2s + 1 ' 

1 s^ + 3s' + 5s + 7 

^a/3W - ^a/3(3) + ^4 ^ + + 2s + 1 " S^ + 2s^ + 25^ + 2s + 1 ' 



2. Incidence graph, pg{l — 1; Z — 1; Z — 1), Z = 4; 5; 7; 8 



N = 2Z^, {6o, &i, &2, &3; ci, C2, C3, C4} = {Z, Z - 1, Z - 1, 1; 1, 1, Z - 1, 1}, 



K — I, K2 — l{l — 1), K3 — l{l — 1), K4 — I — 1. 



Then, 



2/2-1 



R. 



2/2 - 1 - / 2/2-1 



/3(/-l) /2(/-l)' 



/ + 1 2/3 + 1 



a/3(3) — -"-0/3(2) 



/3(/-i) i^i-iy 



1 2/2 + 1 



/3 /2(/-l)' 



3. Hadanicird graph 



A^ = 167, {60, &2, ^3; ci, C2, C3, C4} = {47,47-1,27, 1; 1,27, 47 -1,47}, where, 76 N 
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K = 4'y, «;2 = 2(47-1), /tg = 47, K4 = l. 



Then, 



R. 



I67-I 127-1 

47 + 1 



87-1 



3272(47-1) 47(47-1)' 

647^-47 + 1 ^ ,1 _ 2 

3272(47 - 1) 3272(47 - 1) ' ^^^'^ ^^^'^ ^ 3272 47 - 1 ■ 



4. Distance-reguleir graphs with classical pairameters. Let F denote a distance-regular 
graph with diameter d > 3. We say F has classical parameters {d, q, a, (5) whenever the 
intersection numbers are given by 



[l + a 



i - 1 
1 





d 




i 




i 




1 




1 


){(3-a 


1 



), i = l,2,...,(i. 



), i = 0, l,...,(i- 1, 



(A-i) 



where. 



J 
1 



:= l + g + + ... + ?■ 



For instance, the d-cube is a distance- regular graph with the classical parameters d, q — 1, 
a = and /3 = 1. As an another example of these types of graphs, one can consider the 
graph Bdip"') which is a type of so-called Dual Polar Graphs [?] and is characterized by 
a — 0, P — p", and 



bi^ ^ . '-, i = 0,l,...,d-l 



q-1 

_q' -1 



, i 1; 2, d, 



(A-ii) 



where, p is a prime number and n e N. For example, consider d — A, and q — 2 then we have 
the graph -64(2) with the number of vertices, intersection arrays and valencies 



N = 2295, {bo, 61, 62, 63; ci, C2, C3, C4} = {30, 28, 24, 16; 1, 3, 7, 15}, 
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K = 30, K2 = 280, K3 = 960, K4 = 1024. 

Then, 

2294 p _ p 566 _ 16623 p _ p ^2 _ 16685 

^"^'^^" 34425' - ^«/3'^' + 34425 x 7 " 240975' " ^"^^'' + 240975 " 240975' 

2 16699 



34425 240975 



5. M22 graph 

= 330, {60, bi, 62, &3; ci, C2, C3, C4} = {7, 6, 4, 4; 1, 1, 1, 6}, 
= 7, /s;2 = 42, = 168, ^4 = 112. 

Then, 

R - ?? -R , 161 _ 164 _ 1 _ 1183 

1 _ 113 



6. CoUineEirity graph gen. dodecagon GD{s; 1) 

AT = s^+2{s^+s'^+s^+s^+s)+l, {bo, bi, 62, &3, &4, &5; ci, C2, C3, C4, C5, ce} = {2s, s, s, s, s, s; 1, 1, 1, 1, 1, 2}, 
K — 2s, K2 — 2s^, Ks — 2s^, K4 — 2s^, = 2s^, kq — s^. 

Then, 

+ 2(s^ + + + s + 1) 



-^a/3(3) = -^a/3(2) + 



+ 2(s5 + s4 + g3 + g2 + 5) + 1 ' 

+ 2(^3 + s2 + s + 1) + 3s^ + 4(^3 + ^2 + ,s + 1) 



S6 + 2(6-" 


+ s'l + S'^ + + s) - 


hi + 2(s5 + + s3 + s2 + s) + 1 




+ 2(s2 + s + 1) 


+ 3s^ + 5s^ + 6(s2 + s + 1) 


56 + 2(s5 


+ + + + s) - 


hi 56 + 2(s5 + + s3 + s2 + S) + 1 




s2 + 2(s + l) 


+ 3s^ + 5s'^ + 76-2 + 8{s + 1) 



S6 + 2(s5 + + s3 + s2 + S) + 1 + 2(5^ + g'^ + + + s) + 1 ' 
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s + 2 _ + 3s^ + 5s3 + + 9s + 10) 

Ra/3(^) - RapW + s6 + 2(^5 + + + + s) + 1 ~ + 2{s^ + + + + s) + V 

1 _ + 3s* + 5s^ + 7s^ + 9s + n 

^a/3(6) - ^a/3(5) + ^6 + 2(^5 + ^4 + ^3 + ^2 + 5) + 1 " + 2{s^ + + + + s) + 1 



7. A distance-regulEir graph with 

N = 2[l + l + — —], {bo, 6l,&2,&3,&4;Ci,C2,C3, 04,05} = {1,1-1,1-0,0, 1; l,o,l-o,l-l,l}, 



where I — 7(7^ + 37 + 1) and c = 7(7 + 1), 7 G N. Then, 

l{l - 1) 

K — I, K2 — , K3 — , K4 — I, K5 — 1. 



Then, 



_ c + 2lc + 2l{l - 1) 
~ Z(c + /c + /(/-l))' 

_ c + /c + 2/(/-l) _ 2(c + /) 

^a/3W ^a/3(i) + _ ^^^^ + + _ i)) ^ + + _ ' 

c ^ 2/(/~ l)(/^-c^) + c=^(/ + l) + c/(/-l) 

^a/3(3) - ^a/3(2) + _ _ ^) " - - c)(c + c/ + /(/ - 1)) ' 

^ + 2[/^(2/^-c^-/) + c/(/ + c)] 

^a^W ^a/3(3) + ;(;_i)(c+;c /(/_ l)(/_c)(c + /c + /(/- 1))^ 
^a/3(5) - KpW + + /(/-!)) ~ l{l-l){l-0){0+l0 + l{l-l)) ■ 



8. Biggs-Smith graph 

= 102, {6o, &i, &2, &3, &4, &5, &6; ci, C2, C3, C4, C5, cq, C7} = {3, 2, 2, 2, 1, 1, 1; 1, 1, 1, 1, 1, 1, 3}, 

K — 3, K,2 = 6, H3 = 12, K,4 = 24, K,5 = 24, Kq = 24, = 8. 

Then, 

101 „ ^ 49 150 D D 23 173 

Rapi^) - ^.,«(^) - Kpi^^ + 1^-1^' ^^P^'^ - + 153 - 153' 

10 183 7 190 4 194 

R-a/iW - Ra/3(3) + Y53 ~ 153' ~ ^"'S'*^ + 153-153' - "'^^"^ + 153 ~ 153' 
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1 195 



153 153 



9. Foster graph 

N = 90, {bo, 6i, 62, h: h: h, h, ci, C2, C3, C4, C5, ce, C7, cg} = {3, 2, 2, 2, 2, 1, 1, 1; 1, 1, 1, 1, 2, 2, 2, 3}, 
K = 3, ^2 = 6, = 12, = 24, = 24, kq — 12, = 6, k^ — 2. 

Then, 

89 43 132 4 152 

17 321 11 653 5 663 

- ^a/J(3) + ^ - - + ^ " ^' ^-P''' ' ^-P''' + 270 " 540' 

R -R , 2 _671 _ 1 _675 
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Figure Captions 

Figure-l:(a) Shows the cube or Hamming scheme H{3, 2) with vertex set V — {{ijk) : 
k = 0,l} and relations Ro = {{{ijk), {ijk)) : {ijk) G V}, Ri = {{{ijk), {i'jk)), 
{{ijk), {ifk)), {{ijk), {ijk')) : i ^ i', /, k ^ k'}, R2 = {{{ijk), ii'j'k)), {{ijk), ii'jk')), 
{{ijk),{ij'k')) : i ^ i',j ^ j',k ^ k'} and R, = {{{ijk), {i'j'k')) : i ^ i',j ^ j',k ^ k'} 
respectively. Its non-vanishing intersection numbers are: p5i — 3, Pn = 2, pjj = P21 — 

2, P?2 = P21 = 3> Pi3 = P31 = 1> P22 = 3> P22 = 3> P23 = P32 = 1> P33 = 1- (b)The vertical 
dashed hues denote the four strata of the cube. 

Figure-2: Shows the octahedron or Johnson scheme J(4,2). 

Figure-3: Shows edges through a and /3 in a distance-regular graph. 
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